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Abstract 

We present a classical approach of a mixture of compressible fluids when each constituent has its 
own temperature. The introduction of an average temperature together with the entropy principle 



c/3 

dictates the classical Fick law for diffusion and also novel constitutive equations associated with 
the difference of temperatures between the components. The constitutive equations fit with results 

1/3 _ 

recently obtained through Maxwellian iteration procedure in extended thermodynamics theory of 

"i/3 

£>V multitemperature mixtures. The differences of temperatures between the constituents imply the 

43 ; 

existence of a new dynamical pressure even if the fluids have a zero bulk viscosity. The nonequilib- 
rium dynamical pressure can be measured and may be convenient in several physical situations as 
for example in cosmological circumstances where - as many authors assert - a dynamical pressure 

(N 

C ■ ■ played a major role in the evolution of the early universe. 
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I. MIXTURES IN RATIONAL THERMODYNAMICS 



In the context of rational thermodynamics, the description of a homogeneous mixture of 
n constituents is based on the postulate that each constituent obeys the same balance laws 
that a single fluid does |l[ . The laws express the equations of balance of masses, momenta 
and energies 

r ape 



Ot 



+ div (p Q v Q ) = t c 



< + div ( Pa v Q ® V Q 



dt 



m, 



[a 



1,2,. ..,n) 



(1) 



d{\p a v 2 a + P a e a ) J/1 \ 
— + div | y-p a v a + p a e a y a - t a v a + q 

On the left-hand side, p a is the density, v a is the velocity, e a is the internal energy, q Q is 
the heat flux and t Q is the stress tensor of the constituent a. The stress tensor t a can be 
decomposed into a pressure part — p a I and a viscous part er a as 

t« = -p a l + cr a . 



We consider first only Stokesian fluids (i.e. tr(<r a ) = 0). 

On the right-hand sides r a , m a and e a represent the production terms related to the inter- 
actions between constituents. Due to the total conservation of mass, momentum and energy 
of the mixture, the sum of production terms over all constituents must vanish, 



r a = 0, ^ m Q = 0, y^e Q = 0. 

a=l a=l a=l 

For the sake of simplicity, we ignore in the following the possibility of chemical reactions 
{r a = 0). 
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The mixture quantities p, v, e, t and q are defined as 

n 

p = ^^p a , total mass density, 



a=l 



1 n 



a=l 



mixture velocity, 



1 

£ = £ i + j-J2 p ^ 

P a=l 



internal energy, 



(2) 



t = — pl+er 7 — ^^(p a u a ® Uq,), stress tensor, 



a=l 



q = qj + ^ p a ( e a + — + -u^] u a , flux of internal 

a=l ^ Pa 2 J 



energy, 



where u Q = v a — v is the diffusion velocity of the component a, p = Xla=iPa * s the total 
I n 

pressure, ej = - y p a s a is the total intrinsic internal energy, qj = y^ =1 q Q is the total 

a=l 

intrinsic heat flux and crj = VJ" j cr a is the total intrinsic shear stress. 
We obtain by summation of Eqs. (pQ), 

f <9p 



ft 



+ div (pv) = 0, 



< + div (pv ® v - t) = 0, 



<9(W + p£) , f/1 , 
— h div < I -pt> + pe ] v - tv + q 



(3) 



0. 



which are the conservation laws of mass, momentum and energy of the mixture. They are 
in the same form as for a single fluid. 

In order to compare the balance equations of mixtures and single fluids, we write Eqs. ([T]) 
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in the equivalent form 



( dp 

+ div (pv) = 0, 



dt 



+ div (pv ® v - t) = 0, 



dj^pv_ + pe) + ^ | A 2 + p£ j v _ tv + q ^ = Q) 



(4) 



+ div (Pfe v b) = °> 



(&=!,. ..,n-l) 



<9t 



+ div (p b v b <g> v 6 - tft) = m fe , 



dt I V 2^ 6 ' Vb ~ q& ' = 66 ' 



where the index 6 runs from to n — 1. 



n 

In this multi-temperature model (MT), used in particular in plasma physics [2], we have 
5n independent field variables p a , v a and T a (a = 1, 2, . . . , n), where T a is the temperature 
of constituent a. To close the system (jlj) of the field equations of the mixture process, we 
must write the constitutive equations for the quantities p a ,e a ,q a ,cr a (a = 1,2, ... ,n) and 
m b , e b (b = 1, . . . , n — 1) in terms of the field variables p a , v a and T a (a = 1, 2, . . . , n). 



II. COARSER THEORIES 



Due to the difficulties in measuring the temperature of each component, a common prac- 
tice among engineers and physicists is to consider only one temperature for the mixture. 
When we use a single temperature (ST), Eq. (j3]) 6 disappears and we get a unique global 
conservation of the total energy in the form (jlj) 3 (see for example [3f|). In a recent paper, 
Ruggeri and Simic [4| discussed the mathematical difference between the ST and the MT 
models when the fluid components are Eulerian gases (q a = 0,cr q = 0). They proved that 
the differential system of the ST model is a principal sub-system [5( of the MT model, and 
for large times, MT solutions converge to ST ones. 
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A further step of coarsening theory is the classical approach of mixtures, in which the 
independent field variables are the density, the mixture velocity, the individual temperature 
of the mixture and the concentrations of constituents. 
In that case system (jlj) reduces to the equations 

^ + pdivv = 0, 



where 



dv 

P ^- dlVt = °' 



de 

p— t grad v + div q = 0, 

dt 



P^ + divJ b = 0, (b = 1, • • • ,n - I) 



d _ d 
dt dt dx 



(5) 



represents the material derivative of the mixture motion, 



p 



are the components' concentrations, and 



(6) 



K a=l 



•a H a 



Pa U a = P a (V Q - Vj 



(7) 



,a=l 



are the diffusion fluxes of the components. 

In the classical approach the stress tensor - as in a single fluid - splits up into the pressure 
isotropic part and the viscosity stress tensor a (for Stokesian fluids this is a deviatoric 
tensor) 

t = -pi + a. 

The system (jSj) determines the field variables p, T, v and c& (b = 1, • • • , n — 1). Consequently, 
we need constitutive relations for e, cr, q and J& (b = 1, • • • , n — 1). 

We consider the pressure p(p, T, c&) and the internal energy e(p, T, c&) as given by the equi- 
librium equations of state as they appear in the Gibbs equations for mixture, viz. 

n-l 

V 

ap — 

6=1 



TdS = de — ^- dp — (pi b — pi ) dct, 



(8) 



where /i Q = fj, a (p, T, q,), with a = 1, • • • , n, denote the chemical potentials of the components 
at equilibrium and S is the entropy density of the mixture. 

The entropy balance law is a consequence of equation (jSJ) and system ([5]). By using arguments 
from the thermodynamics of irreversible processes (TIP) presented in js] and 3] chapter 5, 
we obtain the classical constitutive equations of mixtures 



2uT> 



D 



q = L grad (—J + ^L b grad y b ^ 

^ ' 6=1 ^ 

J a = L a grad f — J - ^£ afc grad f 
^ ' b=i ^ 



(9) 



where D D denotes the deviatoric part of the strain velocity tensor D = | (Vv + (Vv) T ) . 
The phenomenological coefficients L, L^, L a and L ab (a,b = 1, • • • , n — 1) are the transport 
coefficients of heat conduction and diffusion. 

Let us note that relation ([9])i is the classical Navier-Stokes equation of a Newtonian (Stoke- 
sian and isotropic) fluid, while 02,3 are generalizations of the original phenomenological 
laws of Fourier and Fick, according to which the heat flux and the diffusion flux depend 
on the gradients of temperature and concentrations respectively (but not on both). The 
TIP permits the temperature gradient to influence the diffusion fluxes and concentration 
gradients to influence the heat flux; both effects are indeed observed and they are called, re- 
spectively, thermo-diffusion and diffusion-thermo or Soret effects. Additionally, the Onsager 
conditions of symmetry yield the following symmetries of coefficients {t| 



J ab 



->ba 



fa, b 



n 



(10) 



and the following inequalities must be satisfied: 



L a L a i, 



and v > 0, 



is a positive definite form, 



(11) 



so that the entropy inequality can be satisfied. 
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III. A CLASSICAL APPROACH TO MULTI-TEMPERATURE MIXTURES AND 
THE AVERAGE TEMPERATURE 



To reveal the relation between the extended and classical models, a formal iterative 
scheme known as Maxwellian iteration is applied (see e.g In the case of the ST model 
the first iterates Jo are calculated from the right-hand sides of the balance laws (jl])5 by 
using "zeroth" iterates - equilibrium values = - on the left-hand sides. The next 

(2) 

step, second iterates 3 a are obtained from the right-hand sides of the same equations by 
putting first iterates on their left-hand sides, an so on. If we apply the first Maxwellian 
iteration the Fick laws of diffusion fluxes ([9]) 3 are obtained. Roughly speaking, the Fick laws 
are obtained by neglecting the accelerations of the relative motions of the constituents and 
the classical theory is an approximation of the ST model (see [3|], Chapter 5). 
In a recent paper, Ruggeri and Simic [8j considered the Maxwellian iteration of system (jl]) 
in the case of a binary mixture of Eulerian fluids; they obtained the Fick laws as a first order 
term of the expansion of the component momentum equations. When each component has 
its own temperature, an additive constitutive equation comes from a limiting case of the 
constituent equation of energy in the form 

6 = L e ( 7l - 72) divv (12) 

where Lq is a new phenomenological positive coefficient and 

6 = T 2 - Ti . 

Eq. (112p is not obtained in classical theory. The aim of this paper is to find a variant form 
of classical approach to recover also equations like (fT2"j) in the general case of mixtures with 
n not necessarily Eulerian compressible fluids. 

In the classical approach, the velocity field v corresponds to an average velocity with 
respect to mass components. Thanks to the Fick Laws, the diffusion fluxes J a = p a (v a — v) 
determine the component velocities v a . 

By analogy with the velocity fields a natural extension of the classical approach is to con- 
sider an average temperature T and Q a = T a — T (a = 1, • • ■ , n) as constitutive quantities; 
similarly with the diffusion velocity fluxes, we name Q a the diffusion temperature fluxes. 
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To define an average temperature T, Ruggeri and Simic assumed that the total intrinsic 
internal energy of the mixture (which coincides with the full internal energy for processes 
not so far from equilibrium (see (j2J)V) is the same in the multi-temperature case as in the ST 
model when the temperature is T [9( . Consequently, T is defined through the local implicit 
solution of the equation 

n n 

P £ = ^2pa£ a {p a , T ) = ^2p a £ a {p a ,T a ). (13) 
a=l a=l 

This choice comes from the case of particular classes of solutions for perfect gases [IJ and 
from the fact that the equation of energy governs the evolution of the common temperature 
T for the ST model. The consequences on the entropy of the mixture will confirm the 
physical grounds of Eq. ffT31 . 

Taking into account Eqs. (j5J) 1; ©4 and (jT31 . Eq. ©3 of energy of the mixture can be 
written as a differential equation for the average temperature T 

dT n ~ l 
P £ ,T—jj- = p 2 £,p div v+ e Ci) div J& + 1 grad v — div q, (14) 

6=1 

where the comma denotes the partial derivative with respect to the subscript. 

As J a (a = 1, • • • , n) are associated with the difference between component and average 

velocities 

<£> a = T a -T (a = 1, • • • ,n) 

corresponding to the difference between component and average temperatures, are non- 
equilibrium thermodynamical variables. Near the equilibrium, Eq. ([TBI can be expanded to 
the first order; then 



where 



cP = ^-( Pa ,T) (16) 

denotes the specific heat at constant volume for the constituent a at equilibrium. Conse- 
quently, Eq. ( TTBT) yields 

n—l 

©n = eoE^ 6 *- 

Pn c V 6=1 
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The definition of the total specific entropy S of the mixture is 

n 

pS = ^2p a S a {p a ,T a ), (18) 

a=l 

where S a (a = 1, • • • , n) are the specific entropies of the components. Let us note that the 
specific entropy S depends only on T and not on 6(, (b = 1, . . . , n — 1). This property comes 
from the Gibbs relation of each constituent, 

T a dS a d&a o"^Pa ' 

pi 

which implies 

OS de 

Ta &F ^ T ^ = &F ( pal ' ^ 



The first-order expansion of Eq. fllSH yields 



P 



Equation ( fT9l) can be evaluated for T a = T; by using Eq. ( TT5l) we deduce 

as. 



a=l 

and consequently 

pS = ^p Q S Q (p Q ,T). (20) 

a=l 

The specific entropy S does not depend on 0& (6 = 1, • • • , n — 1). 

On the contrary, a first-order expansion of the total pressure of the mixture p = 
Y^=iPa{p a ,T a ) together with Eq. ([TTD yields 

p = Po + 7r , (21) 

where 

n n—1 

Po = ^2p a (p a ,T), 7r e = ^V b 6 6 (22) 

a=l 6=1 

and 

r b = y^jnj {p^v^ (Pb,T) ~Pb C V^f (Pn« T )} (& = 1, ■ ■ ■ , n - 1). (23) 

Therefore, the total pressure p of the mixture is a sum of the equilibrium part po depending 
on T and a new dynamical pressure part (as a non-equilibrium term) irg due to the difference 
of temperatures between the constituents. 



We emphasize that the fact that the entropy density depends only on T and not on 
the a justifies - in the present theory - the consideration of O a as constitutive quantities 
(undetermined quantities). This is a key point in the model. 

Consequently, the aim of the next section consists in determining the constitutive equations 
for the a by using the entropy principle and TIP arguments. 



IV. ENTROPY PRINCIPLE AND CONSTITUTIVE EQUATIONS 

We still assume that the internal energy e(p, T, c&) and the equilibrium pressure po(p, T, c&) 
satisfy the Gibbs equation 

n-l 

TdS = de — dp — (p b — p ) deb . (24) 

F 6=1 

The differences between Eq. flE]) and Eq. (124j) consist in the fact that in Eq. (1241 T means 
the average temperature when each component a has its own temperature T a and po takes 
the place of p . At equilibrium, the component chemical potentials are 

Pa (P 5 V) 

Pa = £ *(p a , T ) + ~ T Sa{p a , T). 

Pa 

As in the classical case with a single temperature, the time derivative can be eliminated 
between Eqs (EJi, ^4 and ( 1T4T) . The second-order terms due to the diffusion velocities 
u a (a = 1, • • ■ , n) are neglected; because we consider a Stokesian fluid, the viscous stress 
tensor is deviatoric and we obtain the balance law in the form 



dS 

P 



dt + diV \ f ( q ~ ~ ^ J 

n-l 



q • grad (^j - £ J 6 • grad (^j^j + i tr (J mec , D) , (25) 
where the mechanical flux (for Stokesian fluids) is 

Jmech = & - TT© I. (26) 

Eq. (|25[) can be interpreted as a balance of entropy, if we consider 

l b ~ Pn) ^6 



1 / 

* = T q "5Z ( ^ 

\ 6=1 
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and 

n-l 



£ = q ■ grad (J^j " X> " S rad (^7^) + k ^ ( Jmech D) (2?) 

as the entropy flux and the entropy production respectively. 

We observe that the entropy production is the sum of products of the following quantities: 

thermodynamic fluxes thermodynamic forces 

heat flux q temperature gradient grad (^) , 

diffusion fluxes J& chemical potential gradients grad ( Mi '~ Mn ) , 

mechanical flux J mec h velocity gradient D. 

In accordance with the case of a single temperature model 10. 11| and chapter 5, in 



TIP near equilibrium, the fluxes depend 
general methodology of the TIP 



inearly on the associated forces (see also for the 



13|): 



- For the heat flux and the diffusion fluxes, we obtain the constitutive equations in the 
form of Eqs ©2,3- 

- For Stokesian fluids, the last term of Eq. (127j) corresponding to the mechanical produc- 
tion of entropy can be written in a separated form as 

1 11 

— tr ( J mech D) = — tr (cr D D ) - — n e div v. 

We obtain the constitutive equation of the viscous stress tensor in the form of Eq. Qi and 
the dynamical pressure part due to the difference of temperatures yields 

n-l 

7r tf = 5^r 6 e 6 = -L ff divv, (28) 

6=1 

where L n is a scalar coefficient of proportionality. 

The production of entropy must be non-negative Q, [jjj] and therefore the phenomeno- 
logical coefficients must satisfy the inequalities (llip and 

K > 0. (29) 

Taking into account that the terms given by Eq. ( 1231) depend on (p b , T), from Eq. ( 1281) we 
deduce that the constitutive quantities a (depending a priori on Vv) must be proportional 
to divv: 

O a = £; a divv (a = 1, • • • , n - 1). 
11 



Let 1 1 Man, 1 1 be the matrix such that k a = Ylb=i ^ ab rb > we have 



n-1 



9 a = -J2 M ab ndivv (a = 1,... ,n-l). (30) 



6=1 



Introducing expression ( l30l) into Eq. ( l28l) . we obtain 

n-1 

L w = M ab r a r b > 0, 

a, 6=1 

and assuming the Onsager symmetry property, M ab = M ba (a, b = 1, ■ • ■ , n — 1), we deduce 
that the coefficients M ab are associated with a positive definite quadratic form. 
Finally, the results are the same as in classical theory, but in addition we get new constitutive 
equations (1301) for the difference of temperatures. 

We have considered the simple case of Stokes fluids. If the fluid is non-Stokesian, the 
Navier-Stokes stress tensor of viscosity is 

<x = A (divv)I + 2z/ D D , 

where A is the bulk viscosity. The stress tensor t becomes 

t = -(Po + vr e ) I + a = -p I + 2 v D D , 

with 

P = p + 7l e + 7T ff . 

The nonequilibrium pressure p — po is separated into two different parts. The first one 
7r CT = — A divv is related to the bulk viscosity and the second one Tig = —L n div v is related 
to the multitemperature effects between components. 

For example, we consider the case of mixtures of perfect gases 

P a = — p a T a , e a = c { y ] T a with c ( y ) = — , 

m a m a (7 Q - 1) 

where k and m a are respectively the Boltzmann constant and the atomic mass of constituent 
a. From Eq. f[2"3"j) . we get 

n = p b Cy ] (7 6 - 7„) 

and consequently 

n-1 

6 a =-^M afePb cf ( 7b -7n) divv. (31) 



6=1 
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The dynamical pressure is 



n-l 



Tie = -L v div v with = ^ M ab p a p b c { y ] c ( y ] (j a - j n ) (j b - j n ) . (32) 

a,b=l 

Let us note that Eq. (I3TI) yields the result of Ruggeri and Simic for binary mixtures: 

_v 

( 

P2 c V 



B = Le (71 - 7 2 ) div v with L e = M n [p x $ + p 2 c<?J . (33) 



V. KINETIC VALUES OF THE PHENOMENOLOGICAL COEFFICIENTS 

While the TIP yields only inequalities for the phenomenological coefficients, they can 
be determined by the fact that our approach is a limiting case of system (jlj) which is in 
agreement with kinetic theory. 

For high temperature plasma physics [2], the phenomenological coefficients appear through 
the production terms m b and e& of Eqs. (j3J) and have explicit expressions. Therefore by using 
the Maxwellian iteration, we can evaluate coefficients L ab and M ab . For a binary mixture of 
Eulerian erases the calculus is done in and and we get 

L n = Tj ^T and L e = r e T (34) 

where r j and re are, respectively, the relaxation times of the mechanical diffusion and the 
multi-temperature diffusion 

1 m 1 + m 2 p x p 2 



T.J 



2T' 12 mini2 p 



rg = 1 (m 1 +m 2 ) 2 / 35 n 

3r' 12 m 1 m 2 /0 1 m 2 (7 2 - 1) + P 2 m i(li - 1) 

where r' 12 is the number of collisions per units of volume and time. The coefficient L n of 
the dynamical pressure Tie in Eq. ( 1321) is directly obtained from Eqs. ([33]) 2 , (134]) and ( 1351) 



A;Tpf (mi + m 2 ) 2 



(7i -7 2 ) 2 



3r' 12 mim 2 [p 1 m 2 (7 2 - 1) +p 2 mi( 7l - l)] 2 

This coefficient is null only if the two gases have the same r y l and 7 2 values. 
The ratio of the two relaxation times is 

r e _ 2 (mi + m 2 )p 

tj 3 p 1 m 2 (^ 2 - 1) +p 2 m 1 (7 1 - 1)' 
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This ratio is always greater than 1. In fact if we assume that r y l 7^ 7 2 and denote by 

7max = m ax (71, 7 2 )> we S et ( c = Pi/p) 

2 (mi + m 2 ) 1 2 1 > 

rj 3 cm 2 + (1 - c)mi (7 max - 1) 3 (7 max - 1) ~ 

(the last inequality is because 7 Q < 5/3 for all ideal gases). 

In contrast with the usual approach, in which we consider diffusion and neglect the difference 
of temperatures between constituents, the relaxation time of the multitemperature diffusion 
is always greater than the relaxation time of classical diffusion. 

From Eq. (154"]) . we can see that if tq is small enough, to obtain Lq and L n values close to 
unity, it is necessary to have a high temperature T. Consequently, in the high-temperature 
cases, the multitemperature diffusion terms cannot be neglected. 



VI. CONCLUSION 



In the case of fluid mixtures with different temperatures, a classical framework allows 
one to determine the novel constitutive equation (130]) for the difference of temperatures 
in addition to the classical Navier-Stokes, Fourier and Fick laws. The difference of the 
temperatures produces a new dynamical pressure term which has a physical meaning and 
consequently can be measured. We point out that the dynamical pressure associated with 
multitemperature fluid mixtures exists even if the bulk viscosity is null. For rarefied gases, it 
is well known that the bulk viscosity A is null (Stokes fluids) and consequently no dynamical 
pressure exists in a classical mixture with one temperature; in a relativistic context, this 
pressure exists but remains very small [3]. Due to the nonzero dynamical pressure even for 
Stokes fluids, we conclude that multitemperature mixtures of fluids have a great importance. 
Perhaps such a model may be used to analyze the evolution of the early universe in which 



a dynamical pressure seems essential 16|, ll7j] - Of course, we have to take into account 



a relativistic framework and the chemical reactions. Nevertheless, the preliminary result 
- a new dynamical pressure exists also in the non-relativistic limit and without chemical 
reactions - shows promise and deserves attention for future work. 

We focus attention on the fact that the present results are obtained thanks to an average 
temperature corresponding to the same internal energy as for the single-temperature model. 
It is important to note that the total energy conservation law yields the evolution of T, the 
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entropy principle can be used because the entropy density S depends only on T and the 
Gibbs equation is always defined by using Eq. (T24"l) . 

Finally if we consider Eq. (1141) for the field variables to be spatially homogeneous, T is 
constant as we expect from physical point of view (see ref. 9]). 
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